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The transport equations for a two-dimensional electron gas with spin-orbit interaction are pre- 
sented. The distribution function is a 2 x 2-matrix in the spin space. Particle and energy conservation 
laws determine the expressions for the electric current and the energy flow. The derived transport 
equations are applied to the spin-splitting of a wave packet and to the calculation of the structure 
factor and the dynamic conductivity. 
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Spin injection and coherent control of spins in vari- 
ous nanostructures represent two principal challenges for 
the field of spintronics. Recently, the amount of spin- 
tronics research has grown up extensively with the ul- 
timate goal of applications to the quantum computing 
and information processing 1 . A number of spin-based 
devices have been designed and studie d 2 ! 3 ! 4 ' 5 . Spin ma- 
nipulation in such devices can be achieved by optical^ or 
electrio2*S*2iiii methods or by ferromagnetic gatingii. A 
controlled coupling between spin and orbital degrees of 
freedom is considered to be a particularly promising tool 
of efficient spin manipulation dating back to the seminal 
proposal by Datta and Das^. 

Spin-orbit interaction in two-dimensional electron gas 
(2DEG) confined at GaAs/AlGaAs, GaN/AlGaN or sim- 
ilar heterojunctions arises because of the quantum well 
asymmetry in the perpendicular [z] direction. The re- 
sulting perpendicular electric field leads to the coupling 
of spin to the electron momentum 12 . The strength of 
this coupling can be experimentally tuned by a gate 
voltage! 3 - 1 !. 

Experimental advances in spin manipulation present 
a certain challenge to develop a proper theoretical de- 
scription for various phenomena related to the spin- 
orbit interaction. In particular, modification of univer- 
sal conductance fluctuations and weak localization has 
been studied in quantum dotsi 5 ^ 1 ^* 1 ^. The phenomenon 
of weak localization has been considered in 2DEG as 
welliii-. The Friedel oscillations in the presence of spin- 
orbit interaction 1 ^ and the ac conductivity and the plas- 
mon attenuation^ are calculated. 

The principal goal of the present paper is to derive gen- 
eral transport equations for the spin-dependent distribu- 
tion function of 2DEG including the effects of spin-orbital 
coupling. We assume that the spin-orbit interaction in a 
two-dimensional electron gas has the form, 



a(a x p y - a y p x ) + (3(a x p x - a y p y ), (1) 



where the first term is the Bychkov-Rashba terrain and 
the second term is the linear Dresselhaus (or anisotropy) 
term present in semiconductors with no bulk inversion 
symmetrjiSi. The expression of Eq. Q corresponds 
to the confinement along the (001) growth direction 22 . 
Hereinafter we neglect cubic Dresselhaus terms. The free 
particle Hamiltonian can therefore be written in compact 



notations as, 

r.2 



H = [- fj]a + ai k aip k , a ik 

2m 



P a 

-a -f3 



(2) 



Here, as usual, &o, a x , a y , a z constitute the set of Pauli 
matrices. We use Latin subscripts for spatial coordinates 
and reserve Greek subscripts for the spin indexes. 

Spectral properties. Before deriving the transport 
equation we describe briefly the spectral properties of 
the Hamiltonian QJ. Its diagonalization is straightfor- 
ward and reveals the existence of two spin-split subbands 
in the electron spectrum, 



dp =£, P + A p , e 2p = £ p - A p , (3) 
&> = f^-M, A p = v /p 2 (a 2 + /3 2 ) + 4a/3p^, (4) 



where 

2m 

and p 2 = p\ + p y is the total electron momentum. The 
eigenf unctions corresponding to the eigenstates J2J) are, 



1 



(5) 



with the upper (lower) sign corresponding to the tpi (i/>a) 
state. The phase factor Xp depends on the direction of 
the electron momentum, 



tanxp 



ap x + (3p y 
ap y + (3p x 



(6) 



For the isotropic 2DEG (j3 — 0) the phase \p coincides 
with the angle between the electron momentum and the 
y-axis. 

Further and more convenient description of the spec- 
tral properties can be obtained by considering the spin- 
dependent retarded Green function, defined as usual by, 

iG*p(x,i>) = 0(t - t')((ip a (x)4( x ') + ^COlM*)»- 

Here we have used the shorthand notation for the space 
and time variables, x — (x, t). In a homogeneous sys- 
tem the correlation functions depend on the relative co- 
ordinates x — x 1 only. Using the above expressions (jjj 
|3J) we can write the expression for the retarded Green 



2 



function of free electrons in the momentum representa- 
tion, G R {e,p) = J dtdx G R (x)e let ~ lpx , which after sim- 
ple transformations takes the form, 



G R (e,p) = G*(e,p) = - 



1 q-q + cosxpg^. - sinxpgy 
2 



e - £i P 
1 (T - cos x p (T x + sin Xp<Jy 



e 2r 



(8) 



with the indexes [i = 1, 2 corresponding to the first and 
second terms, respectively, in the last expression of Eq. 
©. 

The central quantity in the transport theory is the den- 
sity matrix, 



f a p{x,x') = ((ip^(x')ip a (x))). 



(9) 



Its value in the thermal equilibrium is related to the 
imaginary part of the retarded Green function via the 
fluctuation-dissipation theorem, 



/(e,p)= ]T n w [G^(e,p)-G R (e,p)}, 



(10) 



here n^p is the Fermi-Dirac distribution function for the 
/i-th state J3J). It is convenient to expand the density 
matrix over the complete set of the Pauli matrices. In 
particular, for the equal-time density matrix, 



/p 



dt 
2^ 



1 



1 



/(e,p) = T^/pO-o + -g P -<3- 



(11) 



we observe according to Eq. (|10fl that in the thermal 
equilibrium, 

/ P = (nip + TJ2 P ), g px = cosx P (nip - n 2p ), ^ 



Spy = - sinx P (wip - n 2f 



Transport equations. In a generic nonequilibrium state, 
the density matrix © obeys a set of conjugated equations 
that can be obtained from the equations of motion for 
the electron operators ip(x) and ij)<{x) determined by the 
Hamiltonian J5J, 



id t 
id t > - 



2m 
V' 2 
2m 



— n + e4 



f(x, x') + iaikViVkfix, x') = 0, 
f(x, x') + ia ik V' k f(x, x')&i = 0. 



(13) 



The equations l|13[) neglect impurity scattering. This is 
justified for ballistic systems when the mean free path ex- 
ceeds the characteristic system size, e.g. in high-mobility 
2DEG in semiconductor heterostructures (we discuss the 
impurity scattering at the end of the paper). In Eq. 
(|13|1 we allowed for the scalar external field 4> x = 0(x, t). 
In the absence of electron-impurity or electron-electron 
collisions no self-energy terms appear in the equations 



(|13fl which makes it sufficient to consider the equal-time 
(t = t') functions only. 

Following the known route of deriving kinetic 
equations 2 ^ we utilize the Wigner transformation for the 
density matrix, 

/ p (x,i) = Jdve-^f(x+ r -,x-^t). (14) 

By taking the sum of equations H13|) in the Wigner rep- 
resentation we obtain, 



[dt + v • V]/ p + iej dq q (/ p _a - / p+ a)e^ x 

+ia i kPk[c'i 1 / P ] + -a ik W k {^u f p } = 0, (15) 

where v = p/m. Here we introduced the spatial Fourier 
transform for the scalar potential <f>q = f dxcf>(x, i)e lqx , 
with the shorthand notation for the momentum integra- 
tion dq = d 2 q/(2ir) 2 . 

Finally, to present Eq. I|15fl in a more transparent way 
we turn to the Pauli matrix representation ljll|l to write, 

[dt + v • V] / p + iej dq - / P +f 

+ a ik V k g pi = 0, (16) 



[dt + v ■ V] g pi + ie J dq q (g p _^ - 5 P +f l )e Jqx 

- [bp x g p ]i + a ik V k f p = 0, (17) 

with the following notation for the precession frequency, 

b pl = 2a ik p k = 2Ap(cosXp,-sinx P ,0). 

Conservation laws. The transport equations (|16I17|) 
are of the Boltzmann type and therefore fulfill certain 
particle and energy conservation conditions which will 
now be obtained. By integrating Eq. IjKil) with respect 
to the momentum we find the continuity equation for the 
particle flow, 



&P+-V-J = 0, 



(18) 



where the electron density and the electric current are 
given respectively by, 

P = J dpf p , j k — e J dp [vkfp + a ik g pi }. (19) 

The terms containing the external potential <^ q cancel as 
is readily seen by the change of integration variables. To 
obtain the energy continuity condition we multiply Eq. 
(|16f) by £ p and Eq. I|17|) by b p and add them together. 
After simple transformations the conservation of energy 
can be written in the conventional form, 



9 t p £ + V-j e =j-E, 



(20) 
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where the energy density and energy current are, 
P e = J dp[£ p f p + b pi g pi ], 



dpvk[£pfp + b pl g pi } + a lk I dp[£ p g pi + b pl f p ] 



Jk 



The equation H20(l means that the local energy change is 
due to the energy flow to the neighboring points in space 
as well as a result of the local Joule heating (right-hand 
side) . 

Wave packet splitting. To give a specific application 
of the derived equations let us now use them to describe 
the propagation of a wave packet in 2DEG with a spin- 
orbit coupling. We neglect a spin-orbit anisotropy (3 = 
for simplicity. The wave packet propagates along the y- 
direction and is uniform along the x-axis. The transport 
equations I|ltjll7fl are then one-dimensional and (with no 
external field applied) yield, 

[d t + vd y ] / = -ad y g x , 
[d t + vdy] g x = -adyf, 
[d t + vdy] g y = -2A p g z 
[d t +vd v ]g z = 2A p g y . (21) 

First, we consider a spin-unpolarized Gaussian wave 
packet injected at the point y — at the time t = 
and moving with the average momentum p, 



/ p (x,i = 0)=a F(y), F{y) = e 



-v 2 s P 2 



5p2 



(22) 



In this geometry the phase factor Xp = 0, which means 
that the precession vector b is directed along the x-axis. 
We also observe that g y = g z =0. The remaining two 
of the equations (|2T)l are easily solved by Fourier trans- 
forming them into a set of linear algebraic equations. The 
general solution of Eqs. (|21|l takes the form, 

f(y,t)=A(y-v + t) + B(y-v-t), 
9x(y,t) = My - v+t) - B(y - v-t), 

where we have introduced subband velocities v± = v±a. 
So far, A(x) and B{x) are two arbitrary functions which 
have to be determined from the initial condition l|21|) 
yielding, A(y) = B(y) — F(y). We find that the incident 
wave packet (|22|) is decomposed into two independent 
constituents oppositely polarized along x-direction and 
moving with different velocities. The spatial distribution 
of the electron density is given by the integral over all 
momenta, i.e. 



p±{y,t) 



1 



2^/ 2 5x{t)^ V[ Sx 2 (t) J ' 

with the average velocities v± — p/m±a, and the Gaus- 
sian width at finite times, 8x 2 (t) = 5p~ 2 + t s \ . To 
observe the spin-orbit induced splitting of a wave packet 
the following conditions should be satisfied, 

dp t 
— < a < — . 

TO dp 



(x - v±tf 



(23) 



The first of the two conditions ensures that the split- 
ting dominates over the wave packet broadening, while 
the second condition means that enough time has to 
elapse before the splitting becomes larger than the in- 
trinsic packet width. 

Now let us consider an injection of a packet initially 
polarized along the y-direction. 



/p(x) = (& + &y)F(y), 



(24) 



The equations for the / and g x components of the density 
matrix remain unchanged with the above analysis still 
valid. The second pair of Eqs. H21fl is independent of the 
first pair and have a solution, 



g y (y,t)=F(y-vt) cos(2A p i), 
g z (y,t) = -F(y - vt) sin(2A p i). 



(25) 



According to the expressions i|25|) the initial spin polar- 
ization precesses with a frequency 2A p around the axis 
perpendicular to the propagation direction. Note that 
the precessing spin propagates with the center-of-mass 
velocity v rather than with the subband velocities v±. 

The above analysis assumes that a wavepacket is in- 
jected with a given momentum p. Such an injection into 
2DEG with a spin-orbit coupling is not easy to achieve. 
For example, injection through an interface with a 'nor- 
mal' (with no spin-orbit interaction) 2DEG2i2£*2& would 
not result in a spatial splitting of a wave packet. This is 
due to the fact that the injection happens with a conser- 
vation of energy rather than momentum. As seen from 
Eqs. ilSI II) the two states with the same energy propagate 
with the same velocity^SiSi within the approximations of 
this paper. However, if we take into account the cubic 
Dresselhaus terms, which have been omitted in our dis- 
cussion, there can be a splitting of velocities at the same 
energy. In order to achieve splitting without the cubic 
terms we need to consider a more complicated setup. As 
a demonstration of principle, we consider the following 
example. Let us inject a wave packet propagating along 
the y-direction with the spin polarized along the interface 
(x-axis), e.g. by injection from a ferromagnetic contact. 
The states forming the wavepacket belong to the subband 
1, with a spin polarization -^(1,1). Let us now switch 
on ac magnetic field along the y-axis rotating the spin 
direction until it is aligned with the z-axis, (1,0), and 
then switch the magnetic field off. The resulting state 
will be an equal mixture of both eigenstates ^=(1,1) and 

^=(1,-1) without any change of momentum (the energy 

is no longer conserved). The velocities of these states are 
different and the packet will split. 

The above picture holds not only for the injection of 
initially polarized packet. If the incident packet is un- 
polarized and has a given energy, upon entering the in- 
terface it will become a mixture of two states: ^=(1,1) 

with the momentum po — ma, and -4^(1,-1) with the 
momentum po + ma. Both velocities remain equal to 
v o — po/m. After switching on the ac magnetic field with 
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the frequency to ~ 2apo (which is a resonant frequency 
for the transition between the two subbands), the first 
state will evolve into the mixture of the states: 1) 
and ^Tg(l) both with the momentumpo + ma, mean- 
ing two different velocities vq and vq — 2a. The same 
reasoning shows that the other initial state will develop 
two velocities vo and vq + 2a. Therefore, the initially 
unpolarized packet will split into three parts. 

Ballistic spin injection. We envisage a spin injection 
from ferromagnetic contacts into ballistic 2DEG among 
the applications for the equations derived above. In this 
case the injection occurs with conservation of energy, and 
can be described by the time-independent solution of the 
equations I|ltjll7fl with the appropriate boundary condi- 
tions, which require a conservation of the normal com- 
ponents of the electric current, Eq. I|19|) . at the inter- 
faces. A corresponding theory would generalize the ex- 
isting approach for the ballistic spin-injection based on 
the ordinary Boltzmann equation 3 -. In the latter case 
the Boltzmann equation method is more convenient for 
the calculation of spin polarization of current and magne- 
toresistance than the direct solution of the single-particle 
Schrodinger equation. 

Structure factor. The electron density fluctuations are 
described by the structure factor— defined as the re- 
tarded correlation function, 

X (x, x') = -i6{t - t'){{p{x)p(x') - p{x')p{x))), (26) 

of the electron density operators p{x) — tpl l (x)'ip a (x). At 
equilibrium the structure factor l|26l) depends on the rel- 
ative coordinates x — x' only. The imaginary part of the 
Fourier transform x(w, q) measures the energy dissipa- 
tion of the external field at a given frequency u> and a 
wavevector q. In the isotropic system the structure fac- 
tor is related to the ac conductivity by the relation, 



a(uj) = lim — — x ( w ,q). 



(27) 



The formula (|27() is readily checked using the Kubo for- 
mula for the conductivity and the continuity equation 

According to the fluctuation-dissipation theorem, the 
structure factor can be determined from calculations of 
the linear response to an external scalar field. The 
field-induced modulation of electron density is related to 
the magnitude of the external perturbation through the 
structure factor according to2&, 



5p(u,q) = ex(w,q)</>(w,q) 



(28) 



The electron density modulation is given by the devia- 
tion of the function / p (i,x) from its equilibrium value, 
5p(u,q) = J dp (5/ p ([j,q), and can be found from the 
linearized equations (|l(jll7fl . In the linear approximation 
by the external field 0(u,q), the distribution function is 
a small deviation 



/p = /p + <5/p> gp = gp + <te P , 



(29) 



from its equilibrium value (|12|) . The linearized transport 
equations (|16I17[) take the form, 

(u - qv)Sf p - a lk q t 8g pk = e0(w,q)(/p + a + / p _|), 
(oj - qv)5g pl - i[h p x (5g p ] 4 - a ik q k Sf p = 



e<Ka;,q)(<? p+?i + g p _ fl , 



(30) 



Solving these equations for the variation of the electron 
density (|19|) we obtain the structure factor with the help 
of the relation (J2EJ , 



lY,f d P [l + (-ir'cos( Xp -x P 0] 



n 



mp- 



Vp+ 



CO ■ 



(31) 



c pp- 



where p± = p ± q/2. The expression (|31[1 with co = 
corresponds to the previously derived result for the static 
dielectric function— To simplify further the subsequent 
discussion we will disregard the anisotropy, = 0, and 
consider the zero-temperature limit T = 0. The two spin- 
orbit subbands are axially symmetric, shown on Fig. 1. 
The subbands are filled up to the same Fermi energy 
level €f but have two different Fermi momenta, p± and 
P2, determined from the equations ei(pi) — 6f, where 
Ci(p) are given by Eqs. II) with (3 = 0. This leads to 
the values, 



Pi = po — ma + 0(m 2 a 2 /pi), 
P2 = Po+ ma + 0(m 2 a 2 /pl), 



(32) 



where po is determined by ep = p 2 ) /2m, namely po is 
the Fermi-momentum in the absence of spin-orbit inter- 
action. Note that the Fermi velocities for the two sub- 
bands, 



dp 



Pa 



ip=Pi 



+ O(m 2 a 2 /p 2 ), (33) 



are the same and (up to higher order terms) equal to 
the Fermi velocity in 2DEG with no spin-orbit coupling 
a = 0. The imaginary part of the structure factor x(cj, q) 
determines the absorption, or Landau damping, of the 
external field at given frequency and wavevector. The 
points in the electron momentum space that contribute 
to the Landau damping correspond to the zeros of the 
denominators. There are total four determined by the 
equations, 



qv ± ap+ ± ap- 



(34) 



with the opposite signs of the last two terms correspond- 
ing to the (gapless) transitions within the same subbands 
[Eq. ©] and equal signs describing the transitions be- 
tween different subbands. 

The terms with p = p' in Eq. I|31|l represent the effect 
of intrasubband transitions. Only indirect (g ^ 0) tran- 
sitions contribute to the imaginary part of the structure 
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Pi P 2 



FIG. 1: Spin-orbit induced subbands of an isotropic two- 
dimensional electron gas, e p = p 2 /2m ± a\p\. At T = all 
states below the Fermi energy are filled. The Fermi mo- 
menta for the two subbands are pi.2 = po T ma. The direct 
transitions, q — (shown by the arrow) are possible for the 
states between the dashed lines, pi < p < P2- 



factor. For small transferred momenta q -C po one can 
disregard the deviation of the cosine factor from unity 
and also approximate rii p _ — rii p+ ~ — qdrii/dp. Taking 
the momentum integral we obtain for the contribution of 
the i-th subband 



UJ 



UJ* 



6(q 



2„.2 



uj 2 ), (35) 



where vq = po/m and vi stands for the density of 
states of the i-th subband at its Fermi surface p = p,: 
vi = &-(l - as). v% = "»(i + as). Note that the 

sum of the two contributions (|35|l is independent of the 
spin-orbit interaction (up to higher-order terms), a con- 
sequence of the fact that the two subbands have the same 
value of the Fermi- velocity. Spin-orbit interaction results 
only in a redistribution of the spectral weight between 
the subbands controlled by the changes in the densities 
of states. 

The terms with \i ^ fj,' in Eq. (|31|l correspond to 
the intersubband transitions. Their contribution to the 
structure factor for ma <C q <C pa is negligible com- 
pared to the above considered intrasubband transitions 
by the factor ~ g 2 /-Po (d ue to the small sin 2 prefactor). 
However, the presence of the two subbands is impor- 
tant as it makes the direct, q = 0, transitions possi- 
ble. The factor n\ v — ni v then defines the momentum 
space available for the direct transitions, p\ < p < P2 
(see Fig. 1), which corresponds to the frequency domain 
2A - 2ma 2 < uj < 2A + 2ma 2 , where A = A p 



J P0 ' 



l'2 



\{uj,q^ 0) 



dp 



aq 

J (w + iO) 2 - 4A 2 , ' 
pi 



(36) 



The imaginary part of this expression is 



S>X(W, 9^0) = 0[Am 2 ^ -(w- 2A ) 2 ]. (37) 



The equation (|36|l corresponds to the previously obtained 
resuh-ii for the optical conductivity <j(uj). The expression 
(|3*f)|l goes to zero with the wavevector, which is easily un- 
derstood by noting that the matrix elements for the tran- 
sitions between ^i(p_) and ip2(p+) states are suppressed 
at small transferred momenta since they are orthogonal 
at q — 0. However, their contribution to the conductiv- 
ity [according to Eq. (|27|l ] remains finite, which is clear 
since the operator of electron velocity has nonzero matrix 
elements for the intersubband transitions even at q = 0. 

The experimental observation of the direct transitions 
(|37|l is feasible in the measurements of the resonant mi- 
crowave absorption in high-mobility semiconductor het- 
erostructures. 

Screened electron- electron interaction and plasmon ex- 
citations. So far our analysis has been restricted to the 
noninteracting electron gas. To incorporate the effects 
of the electron-electron interaction in the random phase 
approximation one has to account for the self-consistent 
electric field induced by the variations of the electron 
density. The potential for this field 4> sc obeys the Pois- 
son equation. In two dimensions the Fourier transform 
of the Poisson equation has the form, 



e0 sc (w,q) = V g p{u),q), 



(38) 



where V q = 2ire 2 /q is the bare Coulomb propagator. The 
random phase approximation (RPA) is then equivalent 
to the substitution (f>(w, q) — > 4> sc {u, q) + q) in the 
right-hand side of Eq. (|30|l . It is straightforward to see 
that the structure factor takes the familiar RPA form, 



xrpa{u, qj 



X(w, q) 



i - Vqx(u, q) 



(39) 



The pole of this expression determines the plasmon spec- 
trum uj — uj q + ij q , where uj 2 = v 2 nq/2, with k = 2-ne 2 v 
standing for the static screening radius. The plasmon 
linewidth is given by the imaginary part of the bare struc- 
ture factor at uj = u>„, 



lq = \VqUq\%X(Uq,q)\- 



(40) 



For the plasmon to be an undamped excitation its fre- 
quency should lie above the electron-hole continuum, 
uj q > qv, which requires k > 2q. As was already pointed 
out in Ref. |2(J the plasmon acquires damping when 
oj q ~ 2 A. Since q -C (qn) 1 ^ 2 ~ ma at this range, the 
direct transitions (|37|l make the principal contribution 
to Eq. l|¥i]l. 

Impurity scattering. The equations presented in this 
paper assume ballistic electron motion. The absence of 
impurities allows one to write kinetic equation as a closed 
set of equations, Eq. JT5J, for the density matrix inte- 
grated over the energy variable e, [see Eq. Ullfl ]. i.e. at 
coinciding times. In the presence of disorder the self- 
energy due to impurity scattering should be added to 
the right-hand side of Eq. I|13|) . In general, since plain 
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waves are no longer eigenstates of the system with im- 
purities, the equations for the distribution function de- 
pending on the momentum p (and not on the energy 
e) become not very convenient. More natural (though 
more complicated) equations would result from integra- 
tion over £ p , similar to the usual spin-degenerate case^i. 
Such equations are beyond the scope of the present pa- 
per. 

Special case a — ±(3. Recently, Schliemann et. ali^I 
proposed a spin field-effect transistor based on a particu- 
lar tuning of the spin-orbit coupling constants such that 
a — (or a = —(3). This special system is expected to 
preserve spin coherence even in the presence of disorder. 
This is due to the fact that the spin eigenstates © are 
independent of the electron momenta, x P = const, there- 
fore a scalar impurity potential does not result in the 
intersubband transitions. The same observation holds 
for the structure factor. Since the matrix elements of 
the density are identically zero for the transitions be- 
tween different subbands the second line of Eq. (|3T)> is 
absent in this case and the structure factor is intact by 
the presence of spin-orbit interaction (up to higher order 



corrections) . 

Conclusions. To summarize, we have derived trans- 
port equations for the distribution function of a two- 
dimensional electron gas with spin-orbit interaction of 
both the Bychkov-Rashba and the Dresselhaus mecha- 
nisms. The distribution function is a 2x2-matrix in the 
spin space. General expressions for the particle and en- 
ergy currents and densities are available in terms of the 
density / p and spin g p distribution functions. The ob- 
tained equations are applied to the wave-packet propa- 
gation in a ballistic 2DEG and to the calculation of the 
density-density correlation function x(w,g). We observe 
that for q > ma the structure factor x(w,g) is almost 
not affected by the spin-orbit interaction, but it reveals 
new features when q <C ma due to the direct transitions 
between different spin-orbit subbands. 
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